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on  the  use  of  simplified  simulation  techniques  and  the  employment  of  generally 
available  components  to  obtain  simulation  of  a  number  of  engine  types  and 
over  a  large  range  of  parameters  at  minimum  cost. 
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Simulator  for  Use  in  Development  of  Jet  Engine  Controls 

Emile  S.  Sherrard 

A  method  of  simulation  and  cost  estimates  of  simulator  components  are  given  for  a 
simulator  to  be  used  in  the  development  of  turbojet  engine  control  systems  for  aircraft. 

The  simulator  employs  typical  d-c  analog  computer  components.  It  is  capable  of  repre¬ 
senting  a  twin-spool,  variable-nozzle,  afterburning  turbojet  engine,  as  well  as  the  engine’s 
controller.  The  simulator  is  of  moderate  accuracy,  flexibility,  and  cost.  It  is  intended 
to  be  used  in  the  determination  of  the  stability  and  performance  of  the  engine  control  system. 

1.  Introduction 


Recommendations  and  cost  estimates  for  a  sim¬ 
ulator  capable  of  representing  a  twin-spool,  after¬ 
burning,  turbojet  engine  and  its  controller  are 
given  in  this  Circular.  The  simulator  is  mtended 
as  a  tool,  to  be  used  by  a  group  of  engineers  engag¬ 
ed  in  developing  engine  control  systems.  The  sim¬ 
ulator  is  designed  to  perform  the  following  tasks : 

(1)  Determine  stability  of  the  engme  control 
system. 

(2)  Determine  performance  of  the  engine  control 
system. 

Its  intended  use  requires  that  this  simulator  be 
of  moderate  accuracy  and  flexibility,  and,  insofar 
as  is  possible,  be  simple  and  low  in  cost.  There 
may,  therefore,  be  large  differences  in  cost,  size, 
and  complexity  between  this  simulator  and  other 
simulators  intended  for;  (a)  Testing  a  complete 
engine  control  [1],^  or  a  component  of  an  engine 
control  system;  (b)  simulating,  testing,  and/or 
developing  a  wide  range  of  engines  and  engine 
control  systems  [2];  or  (c)  use  in  both  engme 
design  calculation  and  control  system  develop¬ 
ment.  [3]. 

For  maximum  usefulness  to  an  engine-controller- 
development  group,  a  simulator  should  have  the 
following  characteristics : 

(1)  The  simulator  must  be  useful  in  both  early 
and  late  stages  of  the  development. 

(2)  The  simulator  must  operate  in  real  time 
in  order  that  both  actual  physical  components 
and  simidated  components  may  be  operated 
together. 

(3)  Transducer  equipment  must  be  provided  so 
that  a  simulated  engine  may  be  operated  either 
with  the  entire  physical  controller,  or  with  some 
of  its  components. 

(4)  Simulator  accuracy  must  be  sufficient  for 
evaluating  control  system  stability  and  perform¬ 
ance. 

(5)  The  cost  of  the  simulator  should  permit  its 
acquisition  by  a  typical  engineering  group  engaged 
in  control  system  development. 

The  first  of  these  requirements  has  been  met  by 
designing  the  simulator  so  that  it  may  easily  be 
expanded  from  a  simplified  representation  of  the 

*  Figures  in  brackets  indicate  the  literature  references  at  the  end  of  this 
Circular. 


control  system  to  progressively  more  detailed 
representations.  This  design  feature  permits  sim¬ 
ulation  with  a  minimum  of  data  in  the  first 
stages  of  the  development  period,  and  full  use  of 
all  available  data  in  the  later  stages  of  develop¬ 
ment. 

Both  analog  computer  and  transducer  com¬ 
ponents  must  have  excellent  dynamic  response  for 
real-time  simulation.  Electronic  computer  com¬ 
ponents  easily  meet  the  dynamic  requirements. 
Servo  mechanisms,  with  an  advertised  frequency 
response  of  40  to  50  cps,  may  be  used  to  obtain 
satisfactory  dynamic  response  for  the  nonelec¬ 
tronic  function  generators  iised  in  simulating  the 
jet  engine.  The  greatest  difficulty  in  obtaining 
adequate  dynamic  response  wdll  be  experienced 
with  the  transducer  ecjuipment.  In  particular, 
the  transducer  for  converting  a  computer  voltage 
representing  speed,  to  the  velocity  of  a  shaft 
capable  of  driving  a  100  hp  load,  will  require  a 
large  improvement  in  its  presently  estimated 
performance  [1].  A  method  for  obtaining  this 
performance  improvement  is  suggested  in  this 
Circular. 

A  great  deal  of  space,  effort,  and  expense  could 
be  saved  by  omitting  transducer  equipment,  and 
relying  on  complete  simulation  of  the  system  by 
analog  computer  components.  This  teclmique  has 
produced  satisfactory  results  for  guided  missile 
systems.  Apparently  it  could  also  be  used  for  jet 
engine  control  systems. 

Accuracy  requirements  for  the  simulation  will 
probably  be  attained  by  the  employment  of  simu¬ 
lator  components  with  an  accuracy  of  1  pei’cent 
or  better.  Such  components  should  give  suffi¬ 
ciently  accurate  representation  of  an  engine  and 
controller  because  design  data  and  experimental 
data  may  have  a  5  to  10  percent  accuracy,  and 
production  tolerances  on  a  particular  engine  or 
controller  may  also  permit  a  5  to  10  percent  varia¬ 
tion  m  important  characteristics.  Analog  com¬ 
puter  components  of  0.1  percent  accuracy  may  be 
used  if  desired,  wuth  onlj’’  a  moderate  increase  in 
simulator  cost.  If  transducer  components  of  this 
accuracy  are  required,  simulator  costs  would  be 
considered  prohibitive  by  most  controller  manu¬ 
facturers. 
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2.  Description  of  Simulator 


Figure  1.  Block  diagram  for  engine-controller  simtdator. 


The  simulator  consists  of  tliree  major  parts,  an 
engine  simulator,  a  controller  simulator  (which 
may  be  replaced  by  a  physical  controller),  and 
transducer  equipment,  to  be  used  when  a  physical 
controller  is  operated  with  a  simulated  engine. 
The  interconnections  between  these  component 
parts  are  slio^wn  in  the  block  diagram  of  figure  1. 
Each  of  these  components  is  described  in  detail 
in  one  of  the  following  sections. 

2.1.  Engine  Simulator 

In  this  portion  of  the  simulator,  typical  analog 
computer  components  are  interconnected  in  such 
a  fashion  that  voltages  in  the  computer  are  made 
to  satisfy  the  equations  that  describe  the  behavior 
of  the  actual  jet  engine.  Inputs  to,  and  outputs 
from  the  engine  simulator  are  in  the  form  of  d-c 
voltages.  Input  voltages,  as  shown  in  figure  1, 
represent:  Main  fuel  flow,  W;  exhaust  nozzle  area, 
N;  and  afterburner  (reheat)  fuel  flow,  R.  Output 
voltages,  as  shovm  hi  figure  1,  represent:  The 
speed  of  rotor  1,  /Si;  the  acceleration  of  rotor  1,  ^i; 
the  speed  of  rotor  2,  S2;  the  acceleration  of  rotor 
2,  A2;  the  turbine  inlet  temperature,  T;  and,  an 
engine  variable,  P,  hereafter  called  engine  pres¬ 
sure,  which  may  be  a  compressor  discharge  pres¬ 
sure,  a  compressor  pressure  ratio,  or  some  other 
hidicator  of  incipient  engine  stall. 

The  engine  outputs  will  respond  to  altitude  and 
aircraft  speed  (ram  pressure)  as  well  as  W,  N,  and 
R.  These  environmental  variables  have  been  held 
constant  in  this  simulator  in  order  to  obtain  sim¬ 
plicity  and  economy  in  engine  representation.  If 
engine  operation,  at  different  values  of  altitude, 
and  ram  pressure  are  desired,  manual  changes  may 
be  made  in  settings  of  the  computer  components 
used  to  represent  the  engine.  Because  there  is  no 
simulation  of  the  d.ynamic  behavior  of  the  airplane 
m  which  the  engine  is  mounted,  provision  for  auto¬ 
matic  changes  of  engine  environment  would  result 
merely  in  convenience  and  rapidity  of  adjustment. 
Because  this  would  entail  greatly  inereased  costs 


and  complexity,  manual  rather  than  automatic 
adjustment  of  engine  environment  is  employed  in 
this  simulator. 

a.  Engine  Equations  in  General  Form 

For  control  purposes,  the  six  equations  given 
below  may  be  employed  to  describe  a  twin-spool, 
variable-nozzle,  afterburning  engine.  The  first 
four  state  that  Ai,  A2,  T,  and  P,  for  constant  alti¬ 
tude  and  ram  pressure,  are  instantaneous  func¬ 
tions  of  five  variables,  the  three  engine  mputs  W, 
N,  and  R,  and  the  instantaneous  rotor  speeds.  Si 
and  S2.  The  last  two  equations  state  that  Ai  is 
the  derivative  of  Si,  and  A2  is  the  derivative  of  S2. 
The  first  four  equations  do  not  involve  any  deriva¬ 
tives  or  integrals  with  respect  to  time.  They  re¬ 
sult  from  the  assumption  that  the  engine  is  a  quasi¬ 
static  system  [4],  i.  e.,  its  thermodynamic  processes 
reach  equilibrium  conditions  in  such  a  short  time 
that  the  fluid  flows,  temperatures,  pressures,  and 
aecelerating  torques  may,  for  control  purposes,  be 
considered  to  have  equilibrium  or  steady-state 
values  at  all  times. 

The  six  engine  equations  are: 


^Ai  (IT,  N,  Si,  S2,  R), 

(1) 

:  =  A2  (IT,  A^,  A'l,  A'2,  R), 

(2) 

’=T(W,N,Si,S2,  R), 

(3) 

=  P  (B^  iV,  R), 

(4) 

(5) 

=  A2. 

(6) 

Computer  mechanization  of  the  last  two  equations 
is  very  simple  on  an  operational  amplifier  type  of 
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analog  computer.  Computer  mechanization  of 
the  first  four  equations  is  difficult  because: 

(1)  Ai,  A2,  T,  and  P  are  not  simple  functions 
of  W,  N,  Si,  S2,  and  R. 

(2)  Complete  data  for  Ai,  A2,  T,  and  P  as  a 
function  of  five  variables  is  usually  unknown  or 
unavailable. 

(3)  There  are  no  commercially  available  func¬ 
tion  generators  that  will  accept  five  different  input 
variables,  and  generate  an  output  that  is  a  com¬ 
plicated  function  of  the  input  variables. 

Mechanization  of  eq  (1)  thi’ough  (4)  on  an 
analog  computer  requires  that  these  equations 
be  replaced  by  a  much  simpler  set  of  equations 
which  still  represent  the  pertinent  behavior  (for 
control  purposes)  of  the  engine.  These  simpler 
equations  must  also  require  a  minimum  amount 
of  engine  data,  and,  if  at  all  possible,  the  data 
required  shoulcl  correspond  to  steady-state  opera¬ 
tion  of  the  engine,  rather  than  to  transient  opera¬ 
tion  of  the  engine.  The  reduction  of  (1)  through 

(4)  to  simpler  equations  satisfying  these  require¬ 
ments  is  discussed  in  the  two  following  sections. 

b.  Method  of  Simplifying  Engine  Equations 

Equations  (1)  through  (4)  are  usually  replaced 
by  equations  that  may  be  mechanized  on  an 
analog  computer  through  the  use  of  either  a 
“component  simulation”  technique  or  a  “trans¬ 
ference  simulation”  technique.  In  component 
simulation,  the  engine  components  (compressor, 
combustor,  turbine,  nozzle,  etc.)  are  individually 
simulated,  and  then  connectech together  to  obtain 
a  simulation  of  the  complete  engine.  This 
technique  reciuires  a  large  number  of  equations 
(two  or  more  for  each  component  simulated),  a 
good  understanding  of  the  thermodynamic  and 
fluid-fiow  characteristics  of  the  simulated  com¬ 
ponents,  and  a  knowledge  of  permissible  approxi¬ 
mations  that  may  be  worthwhile.  This  technique 
simplifies  the  simulation  problem  by  replacing  eq 

(1)  tlu'ough  (4),  which  involve  functions  of  five 
independent  variables,  by  a  larger  number  of 
equations  Avhich  may  involve  functions  of  no 
more  than  two  variables.  References  [2,  3]  give 
examples  of  component  simulation. 

Examples  of  transference  simulation  for  single¬ 
spool  engines  are  fm'nished  by  references  [1,  5]. 
This  technique  simulates  only  the  dynamic  res¬ 
ponse  of  the  engine’s  output  (controlled)  variables 
to  the  engine’s  input  variables  (output  variables 
of  the  controller).  It  does  not  simulate  engine 
variables  that  are  neither  engine  outputs  nor 
engine  inputs.  Transference  simulation  requires 
much  less  background  in  thermod3mamics  and 
engine  design  than  is  required  for  component 
simulation.  It  requires  an  excellent  knoAvledge 
of  useful  approximations  for  describing  the  con¬ 
trolled  responses  of  the  engine.  The  equations 
used  are  fewer  in  number  than  those  employed 
for  component  simulation,  but  for  complex 
engines  they  may  involve  functions  of  more  than 
two  variables. 


Transference  simulation  has  been  employed  in 
this  simulator  for  the  following  reasons; 

(1)  A  transference  simulator  would  be  pre¬ 
ferred  by  users  vdth  a  control-system  background. 
Such  users  would  find  it  easier  to  understand, 
modify,  or  check.  It  would  also  avoid  the  need 
for  intensive  study  of  engine  component  design 
or  operation  by  persons  principally  concerned 
with  engine  control. 

(2)  A  transference  shnulator  will  require  less 
data  than  a  com.ponent  simulator.  This  maj''  be 
only  a  small  advantage  if  the  required  data  is  ob¬ 
tained  from  calculations  made  with  a  digital  com¬ 
puter.  If  the  required  data  is  obtained  from  ex¬ 
perimental  testing,  this  advantage  is  significant. 

(3)  A  transference  simulator  for  engine-con¬ 
trols  development  is  expected  to  cost  less  and  re¬ 
quire  fewer  components  than  a  component  simu¬ 
lator.  This  statement  is  certainly  true  when  a 
sunple  engine  is  to  be  simulated.  It  is  believed 
also  true  for  the  relatively  complex  twin-spool, 
afterbrnning  engine  considered  in  this  Circular. 

c.  Simplified  Engine  Equations 

The  first  sunplified  foim  for  eq  (1)  through  (4) 
is  obtained  by  linearizing  them,  about  points  of 
steadj^-state  operation  (equilibrium  points),  i.  e., 
operating  conditions  of  the  engine  in  which  all 
inputs  are  constant,  rotor  speeds  are  constant,  and 
rotor  accelerations  are  zero.  This  first  step  re¬ 
places  (1)  through  (4)  b}^  four  simpler  equations, 
each  of  which  is  the  sum  of  five  terms.  These 
five  terms  are  then  reduced  to  two  by  a  proper 
choice  of  the  equilibrium  points,  as  detailed  in 
appendix  A.  The  sunplified  equations  for  the 
engine  may  be  written: 


(7) 

(8) 

(9) 

(10) 

Jl 

'AS  1^ 

(5) 

_  A 

dt 

(6) 

above  equations,  5IT=IT— 1 

IT'  and 

882=82—82,  where  IT  and  82  represent,  respec¬ 
tively,  the  instantaneous  values  of  m.ain  fuel  flow 
and  speed  of  the  second  rotor,  IT^'  written  as 
IT' (A’,  ^1,  R)  and  82  written  as  S2{N,  Si,  R)  are 
the  steady-state  fuel  flow  and  steady-state  speed 
corresponding  to  equilibrium  operation  with 
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nozzle  area  N,  speed  of  the  first  rotor  Si,  and 
afterburner  fuel  flow,  B.  All  the  partial  deriva¬ 
tives  are  evaluated  at  this  same  equilibrium  point 
denoted  by  {N,  Si,  R)  and  their  values  will  vary 
as  the  point  {N,  Si,  R)  varies.  In  (9)  and  (10), 
T'=T'iN,  Si,  R),  and  P'=P'{N,  Si,  R)  are 
used  to  denote  the  equilibrium  temperature  and 
equilibrium  pressure  corresponding  to  (N,  Si,  R). 
Equations  (5)  and  (6)  are  repeated  here  for  the 
sake  of  completeness. 

d.  Usefulness  of  Simplified  Engine  Equations 

Equations  (5)  through  (10)  accurately  describe 
the  engine  for  small  values  of  8W  and  882-  Thus, 
they  m.ay  be  employed  in  conjunction  with  cor¬ 
responding  controller  equations  to  determine 
stability  when  engine  and  controller  variables 
execute  only  small  departures  from  any  steady- 
state  operating  point.  Experience  indicates  that 
the  engine  control  system  is  stable,  if  it  is  stable 
for  such  small-departure  operation.  Therefore, 
these  small-departure  equations  provide  an  ade¬ 
quate  description  of  the  engine  for  all  stability 
studies  of  the  engine  control  system. 

These  equations  are  also  useful  in  predicting  the 
dynamic  performance  of  the  engine-controller 
system  for  small  or  slowly  varying  pilot  com¬ 
mands. 

Throttle- burst  operation  is  the  principal  large- 
departure  performance  requiring  simulation.  Such 
operation  occurs  whenever  a  large  and  sudden 
advance  of  the  engine’s  throttle  is  made.  For 
such  advances  (e.  g.,  a  change  from  idle  speed  to 
top  speed)  Ai,  A2,  T—T',  and  P—P'  will  vary 
with  8W  and  5*S'2  in  a  nonlinear  manner  rather 
than  the  linear  manner  indicated  by  eq  (7) 
through  (10).  These  linear  equations  require 
correction  if  the  simulator  is  to  yield,  in  response 
to  a  throttle  burst,  the  maxim.um  engine  acceler¬ 
ations,  maximum  engine  temperature,  maximum 
engine  pressure,  and  time  for  the  engine  to  attain 
a  new  steady-state  operating  condition.  For  a 
particular  equation,  the  amount  of  correction 
required  depends  on  the  magnitude  of  511^  and 
8S2.  The  correction  method  used  in  this  sim.u- 
lator  replaces  6IF,  in  (7)  through  (10)  by 
8W—Ki{8Wy  {i=l,  2,  T,  P).  Ki  is  chosen  so  as 
to  make  the  particular  equation  correct  for  the 
largest  value  of  8W  permitted  by  the  engine’s 
acceleration  boundarj^,  which  is  determined  by 
compressor  stall,  or  maximum  temperature,  or 
rich  blowout.  No  correction  is  made  for  large 
values  of  8S2,  because  it  is  expected  that,  for  most 
engines,  8S2  will  remain  small  during  throttle- 
bm-st  operation.  Details  of  the  8W  correction 
procedure  are  discussed  in  appendix  A.  When 
corrected  in  this  fashion,  the  engine  equations 
may  be  used  to  determine  performance  during 
throttle-bmst  operation. 

e.  Mechanization  of  Engine  Equations 

Equations  (5)  and  (6)  can  be  simply  mecha¬ 
nized  on  an  analog  computer.  Integrating  ampli¬ 
fiers  whose  respective  inputs  are  Ai  and  A2  are 


employed  to  yield  Si  and  S2,  respectively,  at  their 
outputs. 

Equations  (7)  and  (8)  will  be  mechanized  by 
summing  two  terms,  each  of  which  is  the  product 
of  a  partial  derivative  and  an  increment  of  either 
fuel  flow  or  speed.  Each  such  product  is  obtained 
as  the  output  of  a  multiplier.  Because  the  partial 
derivatives  are  functions  of  the  three  variables 
N,  Si,  and  R,  they  may  be  generated  by  a  three- 
variable  function  generator  (3VFG)  whose  inputs  I 
are  N,  Si,  and  R.  The  8W  and  8S2  increments  , 
may  be  generated  by  summing  W  and  —  W'  and  | 
S2  and  —S2,  respectively.  W  is  obtainable  as  1 
an  input  to  the  engine  simulator,  S2  as  an  output 
of  the  simulated  engine,  —  W'  as  the  output  of  a 
3VFG  with  inputs  N,  Si,  and  R,  and  —S2  as  the 
output  of  a  3VFG  with  inputs  N,  Si,  and  R.  The 
block  diagram  for  mechanizing  (7)  and  (8)  is 
given  in  figure  2.  Also  shown  are  two  integrators 
whose  outputs  are  Si  and  S2,  and  whose  inputs 
are  Ai  and  A2.  These  integrators  mechanize 
eq  (5)  and  (6).  Further  details  of  the  mechaniza¬ 
tion  of  eq  (5),  (6),  (7),  and  (8)  for  three  special 
operating  conditions  are  contained  in  appendix  A. 

The  mechanization  of  temperature,  eq  (9),  is 
the  same  as  for  (7)  and  (8)  except  that,  in  addi¬ 
tion,  T' ,  generated  as  the  output  of  a  3VFG  with 
N,  Si,  and  R  inputs,  is  added  to  the  terms  com¬ 
posed  of  a  product  of  a  partial  derivative  and  an 
increment.  The  block  diagram  for  mechanizing  i 
the  temperature  equation  is  given  in  figure  3. 
Further  details  are  contained  in  appendix  B. 

The  mechanization  of  pressure,  eq  (10),  involves 
exactly  the  same  details  as  the  mechanization  of 
eq  (9).  Figure  4  gives  the  block  diagram  for 
mechanizing  the  pressure  equation. 

Figures  2,  3,  and  4  are  consolidated  in  figure  5, 
which  shows  the  block  diagram  for  the  entire 
engine  simulator.  Table  1  lists  required  compo¬ 
nents  and  estimated  costs. 

Correction  of  eq  (7)  through  (10)  for  throttle- 
burst  operation  is  illustrated  in  figure  6  for  eq  (7). 
The  figure  shows  that  (c)Ai/dIF)  8W  has  been 
replaced  by  (bAi/dW)  8W-Ki{8W)\,  where  -Ki 
is  obtained  as  the  output  of  a  function  generator 
whose  inputs  are  N,  Si,  and  R.  During  throttle- 
burst  operation,  the  characteristics  of  some  en¬ 
gines  and  engine  controllers  may  cause  N  or  R,  or 
both,  to  be  constant  or  zero.  Again  —Ki  may  be 
insensitive  to  one  or  both  of  these  variables,  and 
so  may  be  considered  constant  or  zero.  In  such 
cases  —Ki  may  be  generated  more  simply  than 
shown  in  figure  6.  Details  of  the  correction  of  eq 
(8)  through  (10)  are  the  same  as  for  eq  (7). 

2.2.  Data  Required  for  Engine  Simulation 

The  minimum  number  of  engine  environments 
for  which  engine  simulation  might  be  required  are 
considered  to  be  sea  level  with  one  or  two  ram 
pressures,  and  a  representative  altitude  with  one  ; 
or  two  ram  pressures.  At  each  engine  environ¬ 
ment  for  which  system  stability  and  performance 
are  to  be  determined,  the  data  listed  below  will  be  ' 
required. 
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W-Main  Engine 


Note  1: 

For  dry  engine  operation,  function  generator  has  form 
shown  in  Fig.  7,  For  quasi-constant  speed  with  afterburning, 
function  generator  has  form  shown  in  Figures  10  and  11.  For 
variable  speed  with  afterburning,  function  generator  has 
form  shown  in  Fig.  12. 

Figure  2.  Block  diagram  for  simulation  of  engine  speed. 


Table  1.  Estimated  computer  components  required  for  simulating  a  twin-spool  engine 


Type  of  equipment 

Number  of  units  required  for  small  departure  simulation 

Additional  units  re¬ 
quired  for  throttle- 
burst  simulation 

Dry  engine  operation 

After  burning 

Quasi-constant  speed 

Variable  speed 

Summing  amplifiers.  . . 

Integrating  amplifiers . 

Potentiometers  . . . 

30 

2 

50 

3 

12 

36  to  60  (3  to  5  for  each 
tapped  servo  poten¬ 
tiometer). 

8 

30 

2 

50 

3 

12 

108  to  180  (3  times 
number  required  for 
dry  engine). 

20 

30 

2 

50 

3 

12 

8 

36  to  60  (3  to  5  for  each 
tapped  servo  po¬ 
tentiometer). 
$150,000. 

6. 

3  to  12. 

Fast  response  servos .  . 

Tapped  servo  potentiometers _ 

Diode  function  generators . . 

Electronic  multipliers. . . 

Two-variable  function  generators. 

Estimated  cost . .  . 

Number  racks  of  computing 
equipment. 

1. 

1  to  4. 

3  to  5  to  12  to  20  (3  to 

5  for  each  tapped 
servo  potentiome¬ 
ter. 

1  to  4. 

$40,000  to  $70,000 

4 

$100,000  to  $150,000 

6 

$5,000  to  $15,000. 

M  to  1. 
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6W  -  Obtained  froTn  Speed  Simulation  Unit 


Notes: 

1,  For  dry  engine  operation,  function  generator  has  form  shown  in  Fig,  7,  For  quasi¬ 
constant  speed  with  afterburning,  function  generator  has  form  shown  in  Fig,  10  and 
Fig,  11,  For  variable  speed  with  afterburning,  function  generator  has  form  shown 
in  Fig,  1 2. 

2, *  Tt^obtain  simulation  of  throttle  burst  operation,  it  may  be  necessary  to  correct 

— 6W  in  the  manner  shown  in  Fig,  6, 

■yj 

Figure  3.  Block  diagram  for  simulation  of  turbine  inlet  temperature. 


6W  -  Obtained  From  Speed  Simulation  Unit 


1,  For  dry  engine  operation  function  generator  has  form  shown 

in  Fig.  7.  For  quasi-constant  speed  with  afterburning,  function 
generator  has  form  shown  in  Fig.  10  and  11.  For  variable  speed 
with  afterburning,  function  generator  has  form  shown  in  Fig,  12, 

2.  *  To  obtain  simulation  of  throttle  burst  operation,  it  may  be 

necessary  to  correct  6W  in  the  manner  shown  in  Fig,  6, 

^  aw  ® 

Figure  4.  Block  diagram  for  simulation  of  compressor  discharge  pressure. 
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Figure  5.  Schematic  diagram  for  engine  simulator. 


Figure  6.  Correction  of  (d4i/5PF)  for  throttle-burst  operation. 


j  a.  Data  Required  for  Small-Departure  Equations 

(1)  Acceleration  Equations.  Mechanization  of 
the  acceleration  eq  (7)  and  (8)  will  require  the 
following  steady-state  characteristics:  Steady- 
state  fuel  flow,  W' ,  as  a  function  of  nozzle  area, 
N,  speed  of  No.  1  rotor.  Si,  and  afterburner  fuel 
flow,  R;  steady-state  speed  of  No.  2  rotor,  S2  , 
as  a  function  of  N,  Si,  and  R. 

The  following  partial  derivatives  will  be  re¬ 
quired: 

dhli  5^2  dyli  ,  d^2 

dw’ m’  W2’ 

These  derivatives  must  be  known  at  the  steady- 
state  operating  points,  as  functions  of  the  three 
independent  variables,  N,  Si,  and  R. 

(2)  Turbine  Temperature  Equation.  At  and 
near  rated  engine  speed,  the  engine  temperature, 
eq  (9),  will  appear  among  the  equations  of  the 


control  system.  For  this  mode  of  operation,  eq 
(9)  will  be  linearized  about  steady-state  values, 
and  the  following  data  will  be  required:  Steady- 
state  turbine  inlet  temperature,  T' ,  as  a  function 
of  nozzle  area,  N,  speed  of  No.  1  rotor.  Si,  and 
afterburner  fuel  flow,  R.  Also,  the  partial  deriv¬ 
atives: 

bT  ,  bT 
and  ^  I 


must  be  known  at  the  steady-state  operating 
points,  as  functions  of  the  three  independent  var¬ 
iables,  N,  Si,  and  R. 

(3)  Pressure  Equation  Requirements .  If  an  en¬ 
gine  pressure  is  used  in  the  steady-state  control 
of  the  engine,  the  steady-state  pressm’e,  P' ,  as  a 
function  of  N,  Si,  and  R  will  be  required.  In 
addition  bPjbW  and  bPlbS2  evaluated  at  steady- 
state  operating  points  will  be  required. 
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b.  Data  Required  for  Throttle-Burst  Simulation 

(1)  Acceleration  Equations.  In  order  that  the 
small-departure  equations  may  be  corrected  for 
throttle- burst  operation,  both  Ki  and  K2  must  be 
furnished  as  a  function  of  N,  Si,  and  R. 

(2)  Turbine  Temperature  Equation.  During 
throttle-burst  operation,  the  engine  controller 
may  act  to  hold  temperature  at  its  maximum  per- 
missable  value.  Because  tight  control  is  required 
to  prevent  excursions  beyond  the  temperature- 
fixed  acceleration  boundary,  it  would  seem  advan¬ 
tageous  to  linearize  temperature  eq  (9)  about  its 
maximum  permissible  value.  In  such  a  case, 
dT/dB  and  bTldS2  would  have  to  be  evaluated 
at  olf-steady-state  operating  points.  Experimen¬ 
tally,  such  evaluation  is  impractical.  It  may  also 
be  expected  that  calculated  values  of  these  deriv¬ 
atives  at  fractional  engine  speeds  and  off-steady- 
state  operating  points  wdl  be  more  inaccurate  than 
similar  values  for  steady-state  operating  points. 

These  difficulties  render  it  advisable  to  linearize 
tlie  temperature  equations  about  a  steady-state 
value.  If  necessary,  this  equation  may  be  cor¬ 
rected  for  large  values  of  bW,  as  were  the  accelera¬ 
tion  equations.  Such  correction  will  require  that 
a  correction  factor,  K^,  be  furnished  as  a  function 
of  N,  Si,  and  E. 

(3)  Pressure  Equation.  During  throttle-bm-st 
operation,  engine  acceleration  may  occur  with  the 
controller  acting  to  hold  that  value  or  function  of 
P,  which  indicates  the  stall  boundary  of  the  com¬ 
pressor.  As  in  the  case  of  the  temperature  equa¬ 
tion,  it  is  considered  advisable  to  linearize  the 
pressure  eq  (10)  around  steady-state  values,  and 
correct,  if  necessary,  by  the  method  used  to  correct 
the  acceleration  equations.  This  will  require  that 
another  correction  factor  Kp  be  given  as  a  function 
of  N,  Si,  and  E.  If  the  terms  of  pressure  eq  (10), 
have  not  already  been  evaluated  for  use  in  steady- 
state  control,  throttle-burst  operation  will  require 
that  P' ,  dP/dH,  and  c)P/d*S2  he  furnished  as  a 
function  of  N,  Si,  and  E. 

2.3.  Description  of  Controller  Simulator 

As  in  the  engine  simulator,  typical  d-c  analog 
computer  components  are  proposed  for  simulating 
the  engine’s  controller.  In  the  complete  system, 
as  shown  in  figure  1,  output  voltages  of  this  simu¬ 
lated  controller  serve  as  input  voltages  to  the 
simulated  engine,  and  output  voltages  of  the  sim¬ 
ulated  engine  serve  as  input  voltages  to  the 
simulated  controller. 

No  particular  difficulties  are  expected  in  simu¬ 
lating  the  controller.  Present  simulation  tech¬ 
niques  are  considered  adequate  both  for  its  lurear 
representation  and,  if  necessary,  for  the  represen¬ 
tation  of  its  principal  nonlinearities.  The  esti¬ 
mated  equipment  requirements,  and  costs,  are 
listed  in  table  2,  which  shows  a  two-variable  func¬ 
tion  generator  used  to  generate  the  acceleration 
boundary  of  the  engine,  and  three  single-variable 
function  generators  used  to  generate  scheduled 
values  of  N,  Si,  and  E  SiS  a,  function  of  throttle 


position.  A  very  complex  controller  might  re¬ 
quire  appreciably  more  function  generators  and 
linear  computing  components  than  are  listed  in 
table  2;  on  the  other  hand,  a  simple  controller 
will  require  appreciably  fewer.  Therefore,  it 
would  be  advisable  to  reestimate  equipment  and 
costs  for  controller  simulation  whenever  a  very 
complex  or  a  simple  controller  is  to  be  developed. 


Table  2.  Estimated  requirements  for  simulating  the  con¬ 
troller  of  an  engine-controller  system 


Type  of  component 

Number 

required 

rinteerators . 

12 

Computing  ampliflers|g^jjjS^jg^,,  _ 

24 

One- variable  function  generators _  _  _  . .  _  _ .  . 

3 

Potentiometers  . .  .  - 

54 

Two-variable  function  generator . . . 

1 

The  above  equipment,  including  power  supplies, 
can  probably  be  obtained  in  two  racks  of  com¬ 
puting  equipment.  For  0.1  percent  accuracy 
components,  estimated  total  cost  is  $30,000.  For 
1  percent  accuracy  components,  estimated  total 
is  $15,000. 

In  the  early  stages  of  development,  simulation 
of  a  linearized  controller  should  be  adequate.  In 
an  intermediate  stage,  a  control  designer  may  wish 
to  deliberately  introduce  a  nonlinear  element  into 
his  speed  controller,  or  his  temperature  controller, 
or  both.  Because,  for  reasons  of  simplicity,  this 
nonlinearity  will  probably  be  a  function  of  only 
one  variable,  two,  or  perhaps  three  function 
generators  for  a  single  variable  will  be  the  princi¬ 
pal  nonlinear  components  required  for  intermedi¬ 
ate  development.  For  the  later  stages  of  develop¬ 
ment,  transducer  equipment  may  be  required  for 
operating  a  physical  controller  with  a  simulated 
engine,  but  may  be  omitted  if  complete  simulation 
of  the  engine-controller  system  is  employed.  If 
reliance  is  placed  on  complete  simulation,  the 
later  development  period  will  require  simulation 
of  the  essential  nonlinearities  of  controller  com¬ 
ponents,  and  this  will  call  for  additional  computer 
equipment.  The  variety  of  presently  used  con¬ 
troller  components  and  the  differing  types  of 
nonlinearities  exhibited  by  mechanical,  electrical, 
electronic,  and  hydraulic  controllers  prevent 
description  of  a  general  simulator  for  nonlinear 
engine  controllers.  Therefore,  no  estimate  of 
equipment  required  for  nonlinear  controller  simu¬ 
lation  can  be  given  in  this  Circular. 

2.4.  Data  Required  for  Controller  Simulator 

a.  Controller  Component  Data 

Data  describing  the  static  and  dynamic  beha,vior 
of  each  controller  component  will  be  required. 
For  initial  development  studies,  this  data  need  not 
be  extensive  and  may  be  of  only  moderate  ac¬ 
curacy.  For  intermediate  development  studies, 
data  will  be  required  that  will  permit  a  computer 
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study  of  the  effect  on  performance  and  stability 
of  such  nonlinear  effects  as  backlash,  deadband, 
variable  gain,  and  saturation.  For  advanced 
development  work,  data  giving  an  accurate 
i  description  of  each  component  over  its  full  range 
of  operation  will  be  required. 

b.  Engine  Data 

Complete  schedules  for  all  variables  scheduled 
;  by  throttle  position  will  be  required.  The  ac- 
■  celeration  boundary  of  the  engine  will  also  be 
required. 

2.5.  Description  of  Transducer  Equipment 

Transducer  equipment  is  required  whenever  a 
i  physical  controller  is  to  be  operated  with  a  shnu- 
lated  engine.  As  shown  in  figure  1,  its  function  is, 

(1)  to  convert  voltage  outputs  of  the  simulated 
engine  to  physical  variables  which  are  used  as 
inputs  to  the  physical  controller,  and  (2)  to  convert 
outputs  of  the  physical  controller  into  voltages 
i  which  are  used  as  inputs  to  the  simulated  engine. 

It  is  not  considered  necessary  or  desirable  to 
transduce  all  of  the  outputs  of  the  simulated  engine, 
and  the  physical  controller.  Ordinarily,  it  is 
convenient  and  economical  to  simulate  thermo¬ 
couple  response  directly.  Afterburner  fuel  flow, 
which  is  very  large,  is  usually  omitted  for  eco¬ 
nomical  reasons.  Other  transducers  are  considered 
both  feasible  and  necessary;  the  minimum  judged 
adequate  for  a  development  simulator  include 
those  for: 

(1)  Converting  the  voltage  representing  engine 
speed  to  the  velocity  of  a  shaft  capable  of  driv¬ 
ing  a  load  of  100  hp  at  rated  engine  speed. 

(2)  Converting  voltage  to  air  pressure. 

(3)  Converting  main  fuel  flow  to  voltage. 

An  ideal  transducer  would  convert  a  voltage  to 
a  physical  variable,  or  a  physical  variable  to  a 
voltage,  instantaneously  and  accurately.  It  should 
also  be  reasonable  in  cost.  Practical  transducers 
exhibit  neither  instantaneous  response,  absolute 
accuracy,  nor  insignificant  cost.  An  acceptable 
transducer  equipment  for  this  simulator  has  been 
developed  by  Vickers,  Inc.,  for  testing  J-34  engine 

3.  Use  of 

As  noted  in  the  introduction,  this  simulator  is 
intended  to  be  useful  throughout  the  entire  con¬ 
troller  development  program.  In  a  typical  pro¬ 
gram,  it  might  first  be  used  for  preparing  proposals 
on  a  new  control  system,  and  shortly  thereafter 
for  investigating  the  relative  merits  of  different 
control  schemes.  As  the  program  proceeds,  the 
simulator  may  be  used  for  determining  system 
stability,  system  performance,  and  optimum,  or 
necessary  values  of  controller  components.  In 
the  last  stages  of  development,  the  simulator 
might  be  used  for  evaluation  of  controller  hard¬ 
ware  prior  to  tests  with  a  prototype  engine. 

This  variety  of  usage  requires  that  the  simula¬ 
tion  be  easily  simplified,  that  it  be  capable  of 


controllers  [1].  Development  has  been  completed 
on  this  unit,  and  the  whole  test  stand,  or  parts  of  it, 
are  available  on  relatively  short-term  delivery. 
This  unit  is  capable  of  producing  the  principal 
static  environmental  conditions  under  which  a 
controller  would  operate.  It  is  also  of  moderate 
size,  moderate  cost,  and  moderate  power  require¬ 
ments. 

That  unit,  if  employed  in  this  simulator,  maj^ 
exhibit  two  deficiencies.  One  is  the  lack  of  a 
voltage-to-air  pressure  transducer,  which  may  be 
remedied  by  the  addition  of  a  recently  developed 
voltage-to-air  pressure  transducer  [2]  whose  pre¬ 
liminary  performance  seems  to  be  excellent.  The 
second,  and  more  serious  deficiency,  is  the  closed- 
loop  frequency  response  of  the  voltage-to-speed 
transducer,  which  is  believed  inadequate  for  this 
simulator.  Design  calculations  indicate  it  will  be 
unsatisfactory  for  input  frequencies  greater  than 
about  1  cps.  If  used  -with  this  simulator,  the 
closed-loop  response  of  tins  transducer  will  require 
considerable  improvement  over  its  calculated 
response,  to  prevent  the  introduction  of  excessive 
phase  lags  into  the  speed  control  system  when  a 
physical  controller  is  operated  with  a  simidated 
engine.  A  satisfactory  frequency  response  ^vill  be 
such  that  it  introduces  no  more  than  a  20°  phase 
lag  at  the  crossover  frequency  of  the  open-loop 
transfer  function  of  the  engine-speed  control 
system.  Crossover  frequencies  for  various  con¬ 
trol  systems  will  vary,  but  it  is  estimated  that  if 
the  crossover  for  the  voltage-to-speed  transducer 
is  made  to  occur  at  10  cps,  instead  of  1  cps  (as 
presently  estimated),  this  transducer  will  prove 
adequate  for  controller  development  purposes. 

Tins  improvement  in  transducer  response  might 
be  accomplished  in  several  waj^s.  One  method 
might  be  to  insert  a  filter,  consisting  of  analog 
computer  components,  between  the  input  to  the 
transducer  and  the  voltage  representing  engine 
speed.  The  frequency  response  of  this  filter  might 
be  chosen  so  as  to  compensate  the  undesirable 
parts  of  the  transducer’s  frequency  response. 
Another  method  is  to  employ  a  model  technique 
[6]  to  obtain  a  satisfactory  transducer  frequency 
response. 

Simulator 

using  “generalized”  engine  data,  and  that  the 
engine  be  simulated  for  several  different  operating 
environments. 

3.1.  Representation  of  Simplified  System 

The  problem  of  simplifying  the  controller  simu¬ 
lation  will  not  be  discussed  here  except  to  say  that 
it  may  be  easily  accomplished  by  techniques 
frequently  employed  by  control  designers.  The 
most  convenient  way  to  simplify  the  engine 
simifiation  is  to  replace  the  3  VFG’s  by  two-vari¬ 
able,  or  one-variable  function  generators.  Thus, 
if  the  effect  of  afterburner  fuel  flow  R  is  slight,  all 
partial  derivatives,  IF',  and  S'2  may  be  considered 
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functions  of  only  two  variables,  N  and  Si.  If,  in 
addition,  the  engine  is  assumed  to  be  always 
operating  near  scheduled  values  of  N  and  Si,  then 
only  their  values  at  scheduled  operating  points 
are  required  for  W',  S'2  and  the  partial  derivatives. 
In  such  a  case,  these  quantities  may  be  regarded 
as  functions  only  of  N  or  only  of  Si.  This  last 
simplification  would  result  in  replacing  each  of  the 
3  VFG’s  by  a  one- variable  function  generator. 
Such  a  simplified  representation  markedly  reduces 
the  requirements  both  for  engine  data  and  for 
computer  equipment,  also  for  time  required  to  set 
up  a  simulated  engine. 

3.2.  Use  of  Generalized  Variables 

If  all  the  physical  variables  in  eq  (5)  through 
(10)  are  suitably  transformed  [7],  they  will  appear 
in  the  equations  in  a  form  called  “generalized”. 
The  mechanization  of  these  generalized  equations 
will  cause  the  simulated  engine  to  operate  in 
generalized  time.  Operation  in  real  time  may  be 
obtained  by  converting  the  engine  accelerations 
Ai  and  A2  to  their  real-time  values,  and  then 
integrating  to  obtain  Si  and  S2  in  real  time.  Si 
and  S2  may  then  be  connected  to  either  the  trans¬ 
ducer  equipment  or  the  simulated  controller. 
However,  Si  and  S2  must  be  reconverted  to  gener¬ 
alized  values  before  they  may  be  used  in  eq  (7) 
through  (10),  when  these  employ  generalized 


values.  These  conversions  merely  require  multi¬ 
plying  by  constants  given  in  [7],  and  are  easily 
performed.  Generalized  values  of  P  and  T  must 
similarly  be  multiplied  by  a  constant  wherever 
P  and  T  is  connected  to  a  transducer  or  controller 
element;  real-time  values  of  P  and  T  must  also 
be  midtiplied  by  a  conversion  constant  before  they 
can  be  used  in  (7)  through  (10),  when  these 
employ  generalized  values. 

3.3.  Simulation  for  Different  Engine 
Environments 

The  simulator  may  be  required  to  represent  the 
engine  at  several  altitudes  and  air  speeds  These 
environmental  changes  may  cause  significant 
changes  in  the  numerical  value  of  the  functions 
in  (7)  through  (10),  and  so  require  a  number  of 
settings  to  be  changed  in  the  engine  simulator. 
These  settings  will  be  changed  manually  rather 
than  automatically,  as  the  simulated  engine  is  not 
to  “fly”  in  an  airplane  traveling  at  varying  speeds 
and  altitudes.  The  number  of  manual  changes 
required  will  not  be  so  large  that  the  time  to  make 
them  would  be  excessive.  Because  controller 
development  concentrates  upon  a  relatively  few 
operating  conditions,  these  changes,  although  a 
nuisance,  will  not  be  a  significant  handicap  for 
this  simulator. 


4.  Estimated  Costs  for  Engine-Controller  Simulator 


Table  1  lists  the  components  required  for  the 
engine  simulator,  and  their  estimated  cost.  The 
components  are  electronic  except  for  the  fast- 
response  servomechanisms  used  in  the  hmetion 
generators.  Components  of  1  percent  accuracy 
should  prove  adequate  for  the  electronic  com¬ 
ponents.  Amplifiers  and  multipliers  of  0.1  per¬ 
cent  accuracy  may  be  used,  if  desired,  with  an 
attendant  increase  in  cost;  0.1  percent  diode 
function  generators  may  not  be  obtainable.  The 
cost  of  the  servomechanisms  will  be  largely  deter¬ 
mined  by  their  speed  of  response.  Servos  with 
an  advertised  closed-loop  frequency  response  of 
40  to  60  cps  should  prove  adequate  for  meeting 
the  requirements.  Each  servo  should  be  capable 
of  mounting  3  or  4  of  the  tapped  potentiometers 
used  in  the  3VFG’s.  Static  accuracy  of  1  percent 
or  better  is  acceptable,  and  should  easily  be  ob¬ 
tained. 

Table  2  lists  the  components  required  for  the 
controller  simulator,  and  their  estimated  cost. 
The  equipment  listed  should  prove  adequate  for 
simulating  a  large  majority  of  engine  controllers. 


Equipment  requirements  for  very  simple  control¬ 
lers  or  extremely  complex  controllers  may  differ 
from  those  given  in  table  2. 

The  cost  of  transducer  equipment  is  estimated  at 
$50,000  to  $100,000,  depending  on  the  amount  of 
instrumentation  and  computer  components  that 
are  purchased  with  the  test  stand.  It  is  assumed 
each  stand  will  be  individually  ordered  with  vary¬ 
ing  deviations  from  the  original  J-34  design.  An 
additional  transducer  for  converting  a  voltage 
to  an  air  pressure  may  have  to  be  procured 
separately. 

The  total  cost  of  engine  simulator,  controller 
simulator,  and  transducer  equipment  for  a  twin- 
spool  engine  with  afterburning  is  estimated  at  a 
minimum  of  $165,000  and  a  maximiim  of  $280,000. 
The  minimum  figure  consists  of  $100,000  for  engine 
simulation  equipment,  $15,000  for  controller 
simulation  equipment,  and  $50,000  for  transducer 
equipment.  The  maximum  figure  consists  of 
$150,000  for  engine  simulation  equipment,  $30,000 
for  controller  simulation  equipment,  and  $100,000 
for  transducer  equipment. 
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The  engine  simulator  will  be  furnished  with 
instantaneous  values  of  IT,  N,  and  R.  It  will  be 
required  to  generate  Si  and  S2.  Steady-state 
engine  data  may  be  expected  to  furnish  IT'  (N'l, 
S'l,  R'),  i-e.,  the  steady-state  fuel  flow  required  at 
a  given  nozzle  area  N' ,  reheat  fuel  flow  R' ,  and 
rotor  speed  S'l.  Also,  knowledge  of  S2  {N'l,  S'l, 
R'),  the  steady-state  speed  of  the  second  spool, 
may  be  expected.  In  eq  (14),  if  the  steady-state 
values  N'l,  S'l,  and  R'  vary,  the  partial  deriva¬ 
tives  will,  in  general,  also  vary.  Thus,  the  partial 
derivatives  may  be  regarded  as  functions  of  N' , 
S'l,  and  R' ,  and  it  will  be  assumed  that  these 
functions  are  known. 

The  mathematical  form  of  (14)  may  be  simpli¬ 
fied  by  choosing  the  steady-state  values  N' ,  S'l, 
and  R'  equal  to  the  corresponding  instantaneous 
values  N,  Si,  and  R.  Such  a  choice  makes  bN= 
Q^bSi  =  bR,  and  (14)  becomes: 

Si, 

||!{s,-S;(Af',  Si, /.”)}■  (15) 


Simulation  of  Engine  Acceleration 


For  small  departures  of  IT,  N,  Si,  S2,  and  R  from 
a  set  IT',  N' ,  S'l,  S'2,  R'  of  steady-state  values, 
Ai  may  be  approximated  by 


*^+ll 

(14) 

where 


5II"=IT-IT',  8N=N-N',  8Si=Si-S'i, 
8S2=S2—S'2,  and  8R—R—R', 

and  where  the  partial  derivatives  are  evaluated  at 
the  steady-state  operating  point  IT',  N' ,  S'l,  S'2, 
R'.  A  similar  ecpiation  may  be  written  for  A2. 


Similarly,  A2  may  be  written: 

-S'" '?')}+ 

H  I  S,-Sa)V',  Si,  )?')]■■  (16) 

A. 2.  Computer  Representation  of  Engine 
Acceleration 

The  further  discussion  of  eq  (15)  and  (16)  is 
conveniently  divided  into  two  separate  cases: 
“dry  engine”  operation  in  which  there  is  no  after¬ 
burning,  and  “reheat”  operation  in  which  after¬ 
burning  occurs. 
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a.  Dry  Engine  Operation 

In  this  case,  R=0=R',  and  W'  and  S'z  become 
functions  of  N'  and  S'l  only.  We  therefore  write: 

(17) 

82=82  {N',  8[).  (18) 

The  function  W'  (N' ,  8'i)  may  be  generated  by 
the  function  generator  shown  in  figure  7.  The 
voltage  output  representing  W'  (N' ,  8'i)  appears 

S’=S,  -  Speed  of  Rotor  #1 
1  1 


Tapped 
Servo  Pot 

N’=N  -  Engine's  Nozzle  Area 

Figure  7.  Two-variable  function  generator  for  generating 

W'  {N',  S'). 


at  the  arm  of  a  tapped  potentiometer.  Each  tap 
of  the  potentiometer  is  supplied  by  a  function 
generator  whose  output  voltage  represents  W' 
{N',  8'i)  for  a  particular  value  of  N' ,  such  as 
di,  (82,  •  •  •,  ds-  The  arm  of  the  potentiometer 
moves  in  correspondence  to  the  nozzle  area,  to 
select  the  output  of  one  of  the  function  generators, 
or  to  interpolate  linearly  between  two  of  the 
function-generator  outputs.  The  generation  of 
a  voltage  representing  82  (N' ,  8'i)  is  accomplished 
in  the  same  manner. 

Because  the  partial  derivatives  are  evaluated  at 
a  steady-state  operating  point,  they  also  are  func¬ 
tions  of  N'  and  S'l.  They  may  be  generated  in  the 
same  manner  as  W'  and  82.  In  figure  8,  the  four 
variables  W' ,  82,  bA^j'dW,  and  d.Ai/d(S'2  are  shown 
as  the  output  of  the  corresponding  two-variable 
function  generators  FG-W' ,  FG-82,  FG-AA^IbW, 
and  FG-bA],l'682.  Inputs  to  these  function  gen¬ 
erators  are  N'  and  aS'(.  In  figure  8,  A^  is  obtained 
by  summing  (d^i/dIT)5lI^and  {i>Ail'd8^b82.  Each 
of  these  terms  is  obtained  as  the  output  of  a  multi¬ 
plier  whose  inputs  are  a  partial  derivative  factor 
and  an  increment  factor.  Each  partial  derivative 
is  obtained  as  the  output  of  one  of  the  function 
generators.  The  increments  are  obtained  by  sub¬ 
tracting  from  the  instantaneous  fuel  input  W,  and 
the  instantaneous  speed  82,  the  steady-state  values 
W  and  82.  The  latter  are  in  turn  obtained  from 
two-variable  function  generators  whose  inputs  are 
N'  and  8[. 

The  acceleration  A2  may  be  generated  in  the 
same  manner  as  Ai.  Figure  9  shows  the  genera¬ 
tion  of  both  A2  and  Ai.  In  addition,  it  shows 
the  integration  of  Ai  and  A2  to  yield  81  and  82 
winch  in  turn  are  fed  back  to  assist  in  the  com¬ 
putation  of  Ai  and  A2.  In  figure  9 ,  the  inputs  to  the 
engine  are  W  and  N.  The  outputs  are  shown 
as  81  and  82,  the  instantaneous  engine  speeds. 


Figure  8.  Block  diagram  for  simulation  of  engine  acceleration. 
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Figure  9.  Generation  of  engine  spool  speeds  for  dry  engine  operation. 


b.  Reheat  Operation 

Two  possible  modes  of  reheat  operation  will  be 
separately  considered.  The  first  occurs  when 
afterburner  operation  is  prohibited  below  rated 
engine  speed.  In  this  case,  steady-state  values 
of  afterburner  fuel  flow  and  nozzle  area  are  such 
that  the  steady-state  engine  speed  and  turbine 
inlet  temperature  are  both  maintained  at  their  rated 
values.  The  second  mode  of  operation  permits 
afterburner  operation  at  steady-state  speeds 
below  the  engine’s  rated  speed.  The  first  mode 
will  be  called  “quasi-constant  engine  speed,” 
the  second  mode  “variable  engine  speed.” 

For  both  modes,  we  shall  choose  N'=N,  R'  =  B, 
and  S'i=Si  in  order  to  obtain  the  simplified  eq 
(15).  Since  R  is  no  longer  identically  zero,  the 
steady-state  fuel  flow  W',  the  steady-state  speed 
S2,  and  the  partial  derivatives  d^i/dlF  and 
dAi/dS2  are  now  functions  of  three  variables, 
N',  R' ,  and  S'l.  Because  the  first  mode  is  the 
simpler  to  simulate,  it  is  discussed  first. 

(1)  Quasi-Constant  Engine  Speed.  In  this  case, 
S[=Si  is  nearly  equal  to  rated  engine  speed, 
Sw  Since  the  engine  speed  control  and  the 
inertia  of  the  rotor  act  to  maintain  Si  nearly 
equal  to  iSio  at  aU  times,  only  small  departures  of 
S'l  from  Sio  need  be  considered.  This  restriction 
allows  IT'  {N' ,  S'l,  R')  and  the  corresponding 
values  of  S'2  and  the  partial  derivatives  to  be 
generated  rather  simply  for  large  ranges  of  N'  and 
R' .  The  method  is  the  same  for  each  of  these 
four  functions.  It  is  described  below  for  W'  {N' , 
S'l,  R'),  and  diagramed  in  figures  10  and  11. 

At  some  constant  value  d  of  N' ,  the  quantity 
W  (13,  S'l,  R')  is  a  function  only  of  S'l  and  R' . 
Because  S'l—Sw  is  small,  we  write  as  a  good 
approximation : 


W'i^,  S'l,  R')=^W'(^,Sio,  R')+^^{S'i-Sio)  (19) 
where 

W'  (jS,  Sio,  R')  is  a  function  only  of  R' ,  and 
where  5TF7dS(  is  evaluated  at  N'  =  ^,  S'i=Sio, 
and  at  R'.  Consequently,  this  derivative  is  a 
function  only  of  R' . 

The  function  W'  (/3,  S'l,  R'  =  0)  has  already 
been  generated  for  the  case  of  dry  engine  opera¬ 
tion.  This  function  may  be  written,  for  Si  near 
'S'lo,  as 

IT'(d,  s;,  R'  =  0)  =  W'(^,  Sio,  R'^0)+D(S'i-Sio) 

(20) 

where  D  is  the  partial  derivative  of  W  (/3,  Si, 
R'  =  0)  with  respect  to  S'l  evaluated  at  *S'i  =  5'io, 
and  IT'  (/3,  Sio,  R'  =  0)  is  a  constant.  Adding  and 
subtracting  eq  (20)  from  (19)  yields: 

IT'  (13,  S'l,  R')=W'  (/3,  Si  R'=0)-hW'  (13,  Sio,  R')- 

iT'  (d,  -Sio,  7?'=o)+|||,'-d||.s;--Sio|-  (21) 

Equation  (21)  shows  that  W'(0,  S'l,  R')  maj^  be 
obtained  by  the  addition  of  the  steady-state  fuel 
for  dry  engine  operation,  W'(l3,  S'l,  R'  =  0),  the 
output  of  a  function  generator  of  one  variable, 

W'(^,  Sio,  R')-fV'(/3,  Sio,  R'=0) 

whose  input  is  R',  and  the  output  of  a  multiplier 
whose  inputs  are  (S'l  —  Sio)  and  the  output  of  a 
function  generator,  (dlT'/diSJ  — Z)),  whose  input 
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Figure  10.  Generation  of  W  (iV'  =  /3.  ,  jSJ^iSio,  R'). 


Figure  11. 


Generation  of  W'  (N',  S'l^Sio,  R') — steady-state  fuel  flow  for  afterburning  with  quasi-constant  speed. 


is  R'.  A  schematic  diagram  is  shown  in  figure  10. 
As  shown  in  the  diagram,  the  function 
W'(p,  S'l,  R')  may  be  generated  by  three  one- 
variable  function  generators  and  a  multiplier. 
The  value  of  W'((3,  S'l,  R')  for  a  particular  value 
of  (8  is  connected  to  a  particular  tap  on  the  servo 
pot,  as  shown  in  figure  11.  The  other  taps  cor¬ 
respond  to  other  values  of  d-  The  servo  then 
moves  in  accordance  with  N'  to  select  the  function 
for  which  A''  =  /3i  .  .  .,  di  or  to  interpolate  between 
the  two  values  of  d  nearest  N'. 

Details  of  the  generation  of  S2,  dAi/dTF,  and 
dAi/d*S'2  are  the  same  as  for  W',  whose  generation 
is  discussed  above.  Except  for  these  additional 
complexities  in  the  function  generators,  the  block 
diagram  is  the  same  as  for  dry  engine  operation. 

(2)  Variable  Engine  Speed.  The  only  change 
in  the  block  diagram  occurs  in  the  function  gen¬ 
erators  which  produce  W',  S2,  bAijbW,  and 
dAi/d»S'2.  If  these  fnnctions  must  be  generated 


over  large  ranges  of  N' ,  S'l,  and  R'  a  two- variable 
function  generator  must  be  connected  to  each  tap 
of  the  tapped  potentiometers  driven  by  the  nozzle 
servo.  A  schematic  drawing  is  given  in  figure  12 
for  W'{N' ,  S'l,  R').  A  suitable  function  generator 
might  be  that  produced  by  G.  A.  Philbrick  [8]. 
Or,  the  particular  form  of  the  engine  data  might 
allow  it  to  be  generated  by  a  special  technique 
given  by  D.  A.  Elliot  [9]. 

It  is  questionable  whether  these  four  functions, 
when  each  is  a  function  of  three  variables,  will  be 
furnished  the  engine-control  designer.  It  is  very 
unlikely  that  they  will  be  completely  defined  by 
test  data.  Unless  these  engine  data  are  available, 
there  is  little  use  in  simulating  the  engine  for  this 
mode  of  operation.  The  difficulty  of  obtaining 
sufficient  data  may  cause  this  complex  simulation 
to  be  very  infrequently  employed.  However, 
for  the  sake  of  completeness,  the  case  of  variable 
engine  speed  has  been  discussed  in  this  section. 
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Figure  12.  Three-variable  function  generator  for  generating 

W'  iN',  S',  R'). 

A. 3.  Simulation  of  Throttle  Bursts 

If  Sir  and  8S2  are  both  small,  an  excellent  ap¬ 
proximation  for  a  typical  engine  acceleration  is, 
from  (15)  and  (16),  given  by: 


(22) 


In  conjunction  with  the  other  small-departure 
equations  for  the  engine-controller  system,  these 
acceleration  equations  may  be  used  to  determine 
small-departure  stability  of  the  engine  control 
system,  i.  e.,  the  engine  and  its  controller  execute 
only  small  excursions  from  steady-state  operating- 
points.  Because  experience  indicates  the  control 
system  is  stable  if  it  is  stable  for  small  departures, 
these  engine  equations  are  adequate  for  determin¬ 
ing  control  system  stability.  They  are  also  ade¬ 
quate  for  determining  responses  to  small  or  slow 
changes  in  thi'ottle  position,  because,  for  small  or 
slow  changes  in  throttle  position,  the  departures 
of  the  engine  from  steady  state  are  kept  small  by 
the  engine  control  system. 

If  large  and  sudden  changes  in  throttle  position 
are  also  to  be  simulated,  eq  (22)  requires  alteration. 
This  small-departure  equation,  if  used  to  simulate 
large  throttle  bursts,  will  yield  a  higher  accelera¬ 
tion  than  actually  occurs  in  the  engine.  There¬ 
fore,  the  “response  time”  (the  time  to  accelerate 
the  engine  from  one  steady-state  operating  point 
to  another  steady-state  operating  point)  for  large 
throttle  bursts  will  be  smaller  than  that  of  the 
actual  engine.  To  obtain  accurate  simulation  of 
this  response  time,  the  small-departure  equation 
for  acceleration  will  be  corrected  for  large  values 
of  8W,  the  “excess  fuel”. 

During  the  large-departure  operation,  resulting 
from  throttle  bursts,  8W  is,  for  the  major  portion 


of  the  large-departure  period,  equal  to  bWm,  the 
maximum  excess  fuel  which  may  be  furnished  the 
engine  without  the  occurrence  of  stall,  blowout, 
or  engine  over- temperature.  The  acceleration 
Ijroduced  by  8Wm  is,  according  to  eq  (22), 
(dA/dlF)  8Wm-  The  actual  engine  acceleration. 
Am,  produced  by  bWm  is  shown  in  figure  13  as  the 
ordinate  of  the  solid  line  curve  (acceleration 
characteristic)  corresponding  to  an  abscissa 
8Wm-  The  horizontal  dashed  line  in  this  figure 
represents  (dA/dlT)  8Wm-  For  small  8W,  the 
acceleration  characteristic  and  (dA/dlF)  bW  are 
in  good  agreement  with  each  other.  However, 
as  dW  increases,  the  slope  of  the  acceleration 
characteristic  decreases  from  dH/dll'^  at  5lF=0 
towards  a  much  smaller  value  at  8W=8Wm- 
The  exact  shape  of  the  acceleration  character¬ 
istic  may  vary  significantly  with  the  operating 
point  of  the  engine.  Accurate  calculations  of  the 
exact  acceleration  characteristic  are  ordinarily 
not  available.  Accurate  experimentally  deter¬ 
mined  acceleration  characteristics  also  are  ordi¬ 
narily  not  available.  Hence,  it  seems  advisable 
to  compute  an  approximate  value  of  Am  by 
correcting  (dH/dlF)  bWm  and  to  make  this  correc¬ 
tion  a  flexible  one. 


Figure  13.  Typical  variation  of  engine  acceleration  with 
excess  fuel. 

For  most  engines,  the  slope  of  the  acceleration 
characteristic  at  bWm  is  greater  than,  or  equal  to 
zero.  The  approximation  of  the  acceleration 
characteristic  by  (dH/dlF)  bW — K{bWy  will  have 
slope  bAjciW  at  bW^Q  and  a  slope  greater  than 
zero  at  bW=bWm  if  K  is  positive  and  less  than 
K  (dH/dlF)/  AVm-  Therefore,  for  throttle  bursts 
the  acceleration  characteristic  will  be  approxi¬ 
mated  not  by  {bAjbW)  bW  but  by  (dH/dlF) 
bW-K  {bWy  where  Q<K<]i  dH/dlF/  bWm- 
The  exact  value  of  K  to  be  used  will  be  determined 
by  calculated  data,  experimental  data,  or  the 
experience  and  discretion  of  the  control  system 
designer.  The  value  of  Am  will,  according  to  the 
value  of  K  used,  be  between  (dH/dll^)  bWm  and 
Yi  (dH/dlF)  bWm-  The  correction  of  (dH/dlF)  5lF 
for  throttle  bursts  is  shown  in  figure  6  for  rotor 
No.  1. 
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Throttle  bursts  may  also  cause  the  value  of  §82 
to  be  large.  682  is  equal  to  82—82  (N' ,  8'i,  R'). 
The  rotor  inertias  prevent  81  and  82  from  changing 
rapidly  during  a  throttle  burst.  The  rapidity  of 
nozzle  area  change  is  limited  by  the  size  of  the 
nozzle  actuator,  the  inertia  of  the  nozzle,  and  the 
forces  exerted  by  the  exhaust  gases  upon  the 
nozzle.  Permissible  changes  in  afterburner  fuel 
flow,  R' ,  are  usually  made  small  during  throttle 
bursts  in  order  that  the  permissible  value  of  excess 
fuel,  dWm,  will  not  have  to  be  varied  with  R' . 
Therefore,  thi’ottle  bursts  will  not  immediately 
produce  large  changes  in  8[,  N',  R',  82  {N',  8'i, 
R'),  82,  and  882=82—82  (N',  8'i,  R').  Hence,  at 
the  beginning  of  a  throttle  burst  transient,  882 
will  be  small.  It  will  also  be  small  near  the  end 
of  this  transient.  If  82  responds  much  more 
rapidly  to  the  excess  fuel  than  does  81,  882 
will  probably  be  positive  during  the  throttle 
burst  transient.  If  82  responds  much  more 
slowly  to  the  excess  fuel  than  does  81,  882  will 
probably  be  negative  during  the  throttle  burst 
transient.  If  the  response  of  both  82  and  81  to  the 
excess  fuel  is  approximately  the  same,  882  will 
probably  remain  small  and  may  reverse  sign  dur¬ 
ing  the  throttle  burst  transient.  The  values  5S2 
may  assume  during  a  throttle  burst  seem  to  vary 
widely  according  to  the  particular  engine  that  is 
being  simulated.  It  is  felt  that  for  the  majority  of 
engines,  both  rotors  will  exhibit  approximately  the 
same  response  to  excess  fuel,  and  that  882  will  be 
small  during  a  throttle  burst  transient  as  well  as 
at  the  beginning  and  the  end  of  this  transient. 
Therefore,  in  eq  (22),  (dH/d»S'2)  882  will  be  used  in 
computing  acceleration  for  both  small-departure 
and  tlu’ottle  burst  operation  of  the  engine. 

A. 4.  Experimental  Determination  of  the 
Partial  Derivatives 

The  general  equation  for  small  departure  ac¬ 
celeration  has  the  form: 


value  of  and  steady-state  data  by  eq  (25) 

below. 

To  obtain  this  relation,  we  may  write: 


5  W=  W-  W' = W-  W"  +  W"-  W'  =  W-  ' 
W"+8W', 

882=82  82=82 — 82  "h‘S'2  — 82=82 —  [ 

8’2'+88'2 


(23) 


where  W'  and  82  are  steady-state  values  corre¬ 
sponding  to  the  steady-state  operating  point  N', 
8'i,  R',  and  W"  and  82'  are  steady-state  values 
corresponding  to  another  steady-state  operating 
point  N",  8'i',  R"  near  N' ,  8'i,  R'.  If,  in  (23), 
W  is  chosen  equal  to  W" ,  82  is  chosen  equal  to 
82  ,  A  is  zero  at  the  second  steady-state  operating 
point,  and  (22)  may  be  written: 


(24) 


From  (24),  we  have: 


M 

5^2 


dA 

dW 


8W'I882. 


(25) 


The  value  of  c)A/c)W  is  assumed  to  be  known  at 
N',  8'i,  and  R'  from  experimental  testing,  8W' 
and  8S2  may  be  evaluated  from  steady-state  data 
by  choosing  convenient  steady-state  increments 
for  8N'  =  N"-N',  88[  =  8['-8[,  and  8R' = 
R" —R' .  For  rotor  1,  (25)  becomes: 


8W' 

8)82  dW  882 


(26) 


Convenient  values  to  choose  are  8R' =  Q  and 
8N'  =  Q,  and  5ll^'  and  882  as  the  increments  result¬ 
ing  from  a  small  value  of  88'i. 

For  rotor  2  (25)  becomes: 


.  dH  ,rT ,  dH 


(22) 


dH2 _ c)A2  8W' 

d  82  dlF  882 


where  the  partial  derivatives  are  to  be  determined 
at  a  steady-state  operating  point,  {N' ,  8[,  R'). 
The  value  of  bAjciW  may  be  determined  by  ex¬ 
perimental  testing.  Experimental  determination 
of  dH/d/S'2  may  be  difficult  or  impractical.  If  this 
is  the  case  dH/d*S'2  may  be  determined  from  the 


where  8R'  and  8N'  have  the  same  values  as  are 
used  in  eq  (26).  As  a  check  and  a  guard  against 
errors  in  determining  8W'  and  882,  (26)  and  (27) 
should  be  used  to  determine  values  of  dAi/d82  and 
()A2/8>82  for  8W'  and  882  equal  to  values  obtained 
by  choosing  8R'  =  8)  —  88'x,  and  8N'7^Q,  but  small. 
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7.  Appendix  B.  Simulation  of  Turbine  Inlet  Temperature 


B.l.  Derivation,  Simplification  and  Mecha¬ 
nization  of  Temperature  Equation 

The  turbine  inlet  temperature  T  is  assumed  to 
be  a  quasi-static  variable,  and  is  taken  as  an 
instantaneous  function  of  IT,  N,  Si,  S2,  and  R. 
We  denote  this  relation  by  T=T  (IT,  iV,  Si,  S2, 
R),  and  use  the  primed  quantity,  T' ,  to  indicate 
a  steady-state  temperature  corresponding  to 
steady-state  values  IT',  N',  S{,  and  R'.  Sub¬ 
tracting  T'  from  T  yields; 

II  (S.-SD+H  (R-It')-  (28) 

As  in  the  case  of  the  acceleration  eq  (14), 
N',  S'l,  R'  are  chosen  equal  to  the  corresponding 
instantaneous  values  N,  S',  and  R.  Equation 
(28)  then  reduces  to: 


If  desired,  the  temperature  equation  may  be 
corrected  for  throttle-burst  operation  by  the 
method  used  for  correcting  the  acceleration 
equations.  The  amount  of  correction  required 
should,  for  most  engines,  be  relatively  small. 

B.2.  Experimental  Determination  of  Partial 
Derivatives 

The  two  partial  derivatives  in  (29)  are  to  be 
evaluated  at  steady-state  operating  points,  and 
should  therefore  be  much  easier  to  determine  by 
experimental  testing  than  if  their  value  were 
required  at  transient  operating  points.  Experi¬ 
mental  determination  of  dT/dIT  vdll  probabl}^ 
be  practical.  Experimental  determination  of 
dT’/d»S'2  may  be  difficult  or  impractical,  in  which 
case  this  derivative  may  be  calculated  in  terms 
of  dT/5ITand  steady-state  data  by  eq  (32)  below. 

This  equation  is  developed  by  taking  the 
difference  of  two  steadv-state  temperatures  T" 
(N",  S'l',  R")  and  T'  (iV',  S'l,  R').  The  values 
8N'=N"-N',  8S'i=S'i'-S'i,  and  8R'=R"-R' 
may  be  chosen  small  so  that  we  may  write  as 
a  good  approximation: 


where 

SW^W-W'=W-W'  (S'l,  N',  R'),  dS2=S2- 

S'2  {S'l,  N',  R'),  (30) 

and  the  partial  derivatives  in  (29)  are  evaluated 
at  a  steady-state  operating  point  determined 
N',  S'l,  and  R' . 

The  mechanization  of  eq  (29)  and  (30)  requires 
generation  of  IT'  (ySJ,  N' ,  R'),  S'2  (S'l,  N',  R'), 
and  the  partial  derivatives  that  are  likewise 
functions  of  N' ,  S'l,  and  R' .  The  generation  of 
ST  is  shown  in  schematic  form  in  figure  3.  It  is 
similar  in  detail  to  the  generation  of  an  engine 
acceleration.  To  obtain  T,  the  instantaneous 
temperature,  T'  is  added  to  hT.  T'  is  also  a 
function  of  A',  S'l,  and  R' ,  and  is  generated  by  a 
function  generator  whose  details  are  identical 
wuth  the  function  generator  which  generates 
W'  (N',  Si,  R'). 


§T'=T"~T'=^y{W"~W')+~  (S''-S'2)  = 
dT  ir/i 

By  a  suitable  choice  of  8N',  SS'i,  and  8R',  and  with 
not  more  than  one  of  them  equal  to  zero,  the  steady- 
state  temperatures  may  be  made  the  same.  For 
such  a  choice,  8T'  =  0  and  (31)  becomes: 

<32) 

where  5 IT'  and  8S'2  are  steady-state  increments 
about  IT'  and  S'2  corresponding  to  steady-state 
increments  8N' ,  8S'i,  and  8R'  about  N' ,  S'l,  and  R'. 

In  (32)  dT/dB^is  supposedly  known  at  N' ,  S'l, 
and  R'  from  experimental  values.  A  convenient 
choice  of  8N' ,  8S'i  and  8R'  wluch  makes  8T'  =  0 
will  permit  5 IT'  and  8S'2  to  be  determined  from 
steady-state  (la.ta. 
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